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Abstract. In this not e we stress the necessity of a c areful check of both the theoretical arguments and the numerical 
experiments used bv lLopez-Caniego et al. I IT20051 to support the superior performances of the biparametric scale 
adaptive filter (BSAF) with respect to the classic matched filter (MF) in the detection of sources on a random 
Gaussian background. 
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1. Introduction 

A popular tool for the detection of sources embedded in 
a noise background is the matched filter (MF). This is be- 
cause its optimal properties. In particular, in the case of 
Gaussian noise, it can be shown that through MF it is pos- 
sible to construct uniformly most powerful tests, i.e., tests 
that for a given probability of false alarm Pfa (the prob- 
ability of a false detection) provide the best probability of 
detection Pq (the probability of a correct detection) . This 
is guaranteed by the Neyman-Pearson (NM) theorem. For 
the sake of clarity and completeness, some arguments are 
presented in the following. 

Let's suppose that x = [x(0), x(l), . . . , x(N— 1)] T is an 
observed one-dimensional discrete signal. In its simplest 
form, the detection problem consists in deciding whether 
x is constituted by a pure zero-mean Gaussian noise n 
with covariance matrix C (hypothesis Hq) or it contains 
also the contribution of a deterministic component s (hy- 
pothesis -Hi). Both C and s are supposed known. In other 
words, the source detection problem is equivalent to a de- 
cision problem where one has to distinguish between the 
two hypotheses: 



Under Hq the probability of x is given by 



p(x\H Q ) = p G {n) = p G {x) 
1 



(2 7 r) JV / 2 det 1/2 (C) 
whereas under H\ it is 



exp 



--x T C~ x x 
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, (2) 



p(x\Hx) = p G (n) = p G (x - s) 



1 



(2 7 r) JV / 2 det 1/2 (C) 
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\(x-s) T C-\x-s) 



(3) 



The choice between Hq and H\ requires the definition of 
the so called critical or acceptance region R*. This is set 
of values in that map into the decision H\ or 



R* = {x : decide iJior reject Hq}. 



(4) 



Once i?* has been fixed, the quantities Pu and Pfa are 
given by 

(5) 



and 



x = n; 
x = n 



(1) 



Pn = I p(x\H{)dx, 

JR., 

p(x\Ho)dx. 



Pi 



(6) 
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The Neyman-Pearson (NP) theorem tells us how to choose 
P* in such a way to maximize Pd when we are given the 
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likelihood functions p(x\H ) and p(x\Hi), and a value a 
for P FA . 

Neyman-Pearson theorem: to maximize Pd for a given 
Pfa = ol, decide Hi if the likelihood ratio 



L(x) 



v(x\Hi) 



p(x\H ) 

where the threshold 7 is found from 



> 1, 



Pi 



FA 



{x:L(x)>7} 



p(x\H )dx = a. 



(7) 



(8) 



By inserting Eqs.|J5J-l(BJ in Eq.JTJ, it is possible to see that 
Hi has to be chosen when for the statistic T(x) (called 
NP detector) it is 



T{x) = x T C~ 1 s > i 



with 7' such as 



Pfa = Q 



(9) 



(10) 



Here, Q[x) = 1 — $(x) with the cumulative dis- 

tribution function of the standard Gaussian distribution. 
Equation (|10l) is due to the fact that, being given by a 
linear combination of random Gaussian quantities, T(x) 
is a random Gaussian quantity too. For the same reason 
it happens that 



Pd = Q Q- 1 (Pfa) 



(11) 



with e = s T s the signal energy and a the standard devia- 
tion of the noise n. In Eq.JjJJ, the quantity C _1 s provides 
the classical MF. Therefore, this filter can be expected to 
have optimal properties. An important point to stress is 
that these properties do not degrade very much when the 
amplitude and/or the position of the source are unk nown 
and have to be estimated from the data themselves IlKav I 



In iLopez-Caniesro et al. I ((2005) (hereafter LHBS ) an 
alternative approach is proposed that is claimed to pro- 
vide filters with superior performances with respect to 
MF. Authors work in the context of one-dimensional con- 
tinuous signal and under the hypothesis of sources with 
symmetric profile. Their idea consists in exploiting the 
joint distribution p(v n ,K n ) of the normalized maxima v n 
and of the corresponding normalized curvature K n that, 
for a continuous Gaussian random noise, is given by the 
Rice's formula 



PR.{v ni K n 



cxp 



2pv n n r , 



V2vr(l - p 2 ) L 2 (1 - P 2 ) 

(12) 

Here, v n = £/er 6 (-00, +00), n n = — C/0-2 S (0, +00) 
with £ the value of the maxima of the signal and £" the 



corresponding second derivative, whereas of and P{q) are 
the moment of order 2i 



/>oo 

2 / q 2l P(q)dq, 
Jo 



(13) 



and the power spectrum of the noise process, respectively. 
Of course, for the first normalized derivative of the max- 
ima it is S„ — £' ' ju\ = 0. 

If a source is embedded in the background and v, 5 
and k are observed quantities, the NP approach provides 
the following decision problem 



Ho 



v = u n , 5 = S n = 0, K = n n ; 
v = v n + v s , 5 = S n + 5 S = 0, n = n r _ 



(14) 

where v n , v s are the values of the noise process and of the 
source in correspondence to the position of the observed 
maximum, whereas 5 n , 8 S and k„, k s are the corresponding 
normalized first derivatives and curvatures, respectively. 
In other words, under Hq it is p(v, 0, k\Hq) = pn{y n , k„), 
whereas under Hi it is p{y 1 0, k\Hi) = pf(jy ni S ni K n ). A 
point to stress is that the likelihood p/(? n , S n , K n ) is not 
given by the Rice's formula. This is because, in general, 
there is no guarantee that the position of the observed 
maximum coincides with the position of the peak of the 
source and/or with the position of one of the peaks in the 
noise process. 



2. Some comments on the LHBS approach 

The LHBS approach deserves some comments. Among 
these, three are of particular relevance. The first one is 
that, instead of pf(v n ,5 n ,Hn), under H\ authors use a 
likelihood in the form 



p(v, k\Hi) 



x exp 



iy 2 



■ s ) 2 + (k - n s ) 2 - 2p(v - v s )(n - k s ) 



2(1 



(15) 



where v G (—00, +00) and k € (0, +00). Actually, it is not 
clear from where this equation comes out. In fact, even in 
the unrealistic hypothesis that the position of the peak of 
the source has to coincide with the position of a maximum 
of the noise process, similarly to Eq.JTJl it should be 

p(v,n\Hx) 



a/2tt(1 - p 2 ) 



x exp 



(y 2 - v s f + (k - n s ) 2 - 2p(v - u s )(k - k s ) 



2(1 



(16) 



with v G (— 00, +00) and k G (/c s ,+oo). LHBS justify 
the use of the likelih ood (|15fl only with a reference to 
iBarreiro et al. I ([2003) . However, also in that work no de- 
tails are provided. 
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The second comment regards the fact that LHBS do 
not work on the original signal x, but on a filtered ver- 
sion of it. Th e reason is that this operation modifies i?* 
(see Fig. 5 in iBarreiro et al. ll2003|) . Hence, LHBS try to 
define a filter for which, once fixed Pfa, R* is such as to 
maximize the corresponding Pp. This is accomplished by 
parametrizing some well known filters (e.g., MH and the 
Mexican hat wavelet) and by computing the values of the 
parameters through an optimization procedure. For exam- 
ple, the biparametric scale adaptive filter (BSAF), claimed 
by LHBS to outperform MH, is a modified version of MF 
where two coefficients are assumed to be free parameters. 
However, authors provide only an heuristic explanation of 
the reason why such an approach should be expected to 
work. No rigorous argument is given. Probably, the conjec- 
ture is that their procedure provides a two-steps (parame- 
ter optimization + filtering) approximation of the solution 
of the decision problem based on the likelihood ratio 



L(x, v, k) = — — — > 7, 

p(x,v,K\H Q ) 



(17) 



that (at least in principle) should better exploit the avail- 
able a priori information. The main problem to work di- 
rectly with such likelihood ratio is in the difficulty to obtain 
the joint distributions p(x, v, k\Hq) and p(x, v, k\Hi) that 
can be quite complex. In any case, the usefulness of the 
proposed procedure is not obvious and should be proved 
through rigorous theoretical arguments. The conclusions 
drawn only on the basis of numerical experiments lack 
generality and can be potentially misleading. This is es- 
pecially true for the discrete case where the quantities v 
and k have to be computed by using the data x. 

The last comment regards the procedure followed by 
authors in the numerical experiments. In particular, it is 
questionable that they investigate the performance of their 
detection test only on the sources that, after the filtering, 
show a peak placed at the position of their true maximum. 
This is certainly an unrealistic condition. Moreover, in this 
way there is the risk to favour the selection of the sources 
that are located in correspondence to a maximum of the 
noise, i.e., those whose likelihood p(v, k\Hx) is given by 
Ea. flfip . In the discrete case, it easy to realize that this is 
very probable when, for example, a source is superimposed 
on a white noise background and no filtering operation is 
carried out. On the other side, BSAF is a symmetric filter 
that works principally on spatial scales smaller than a pix- 
els. Hence, the necessity to check for the above mentioned 
risk. The correct procedure requires that the detection test 
be applied also to the sources that after the filtering show 
a maximum at a position different from the true one. 



As last comment, it is useful to stress that the BSAF 
corresponding to a specific problem can be obtained only 
through a very involved procedure based on stringent a 
priori assumptions (e.g., Gaussianity of the background, 
symmetry of the source profile) . All this makes the proce- 
dure proposed by LHBS quite rigid. In this respect, MF 
certainly represents a much more flexible solution. 
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3. Conclusions 

From the discussion presented above, it is evident the ne- 
cessity of a careful check of both the theoretical argu- 
ments and the numerical experiments used by LHBS be- 
fore BSAF can be claimed to be superior to the classic 
MH in source detection problems. 



